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Abstract— Three software systems, ERCS08, ISICS 2011 and 
Smit's code, that implement theoretical calculations of inner shell 
ionization cross sections by proton impact, are validated with 
respect to experimental data. The accuracy of the cross sections 
they generate is quantitatively estimated and inter-compared 
through statistical methods. Updates and extensions of a cross 
section data library relevant to PIXE simulation with Geant4 are 
discussed. 

Index Terms — Monte Carlo, simulation, Geant4, ionization, 
PIXE 



I. Introduction 

THE calculation of inner shell ionization cross sections by 
proton and ion impact is an important component of the 
simulation of FIXE (Particle Induced X-ray Emission) and the 
analysis of experimental PIXE spectra. The ECPSSR (Energy- 
loss Coulomb Perturbed Stationary State Relativistic) [1] the- 
ory with its variants is regarded as the standard approach for 
cross section calculations in the domain of PIXE applications, 
which typically concern the energy range up to a few tens MeV 
and the whole range of elements in the periodic system. It 
provides inner shell ionization cross sections for PIXE analysis 
codes such as GeoPIXE [2], GUPIX [3]-[5], PIXAN [6], 
PIXEF [7], PIXYKLM [8], Sapix [9], and TTPIXAN [10], 
and for specialized PIXE simulation codes [11]-[13]. 

Several cross section models for the computation of inner 
shell ionization by proton and a particle impact are available 
in a package for PIXE simulation [14] released in Geant4 
[15], [16] 9.4 version; they include models based on the 
plane wave Born approximation (PWBA) [17], the ECPSSR 
model in a number of variants and a collection of empirical 
models, deriving from fits to experimental data. The PWBA 
and ECPSSR cross section models (with variants) exploit 
tabulations produced by the ISICS [18] code for K, L and 
M shells, which have been assembled in a data library [19] 
publicly distributed by RSICC (Radiation Safety Information 
Computational Center at the Oak Ridge National Laboratory). 

A new version of ISICS and an entirely new code for 
the calculation of ECPSSR cross sections (with variants), 
ERCS08 [20], have become available since the publication 
of the previously cited paper on Geant4 PIXE simulation. 
This paper evaluates the cross sections deriving from these 
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evolutions, along with those calculated by Smit's code [21] 
(identified in the following as KIO-LIO), which is consistent 
with pristine ECPSSR formulation. 

The K and L shell proton ionization cross sections produced 
by these three theoretical generators are compared with refer- 
ence collections of experimental data to assess their validity, in 
compliance with the IEEE Standard for Software Verification 
and Validation [22]. The results of this validation process 
document quantitatively the relative merits of the three codes, 
evaluate the impact of the newly available calculations on 
Geant4 accuracy and identify the state-of-the-art of theoretical 
cross sections for PIXE simulation with Geant4. 

II. Theoretical overview 

In the PWBA approach [17], the first-order Born approxi- 
mation is used in scattering theory to describe the interaction 
between an incident charged particle and an atomic target. This 
treatment is justified when the atomic number of the projectile 
is much smaller than the atomic number of the target, and the 
velocity of the incident particle is much larger than the velocity 
of the target-atom electron velocities. 

The PWBA cross section in the center of mass system for 
the ionization of a given shell is given by 

(TPWBA = (T0^"'FQ,^^) (1) 

where: ^ 

f^o - Snal (2) 

flQ is the Bohr radius, Zi is the projectile atomic number, Z2 
is the effective atomic number of the target atom, F is the 
reduced universal function, and the reduced atomic electron 
binding energy and reduced projectile energy are given by 
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with E, M and U representing the energy, mass and atomic 
binding energy of the projectile and the target, respectively 
identified by the indices 1 and 2. The analytical formulation 
of the reduced universal function F can be found in [18]. 

The ECPSSR theory [1] was proposed to address the short- 
comings of the PWBA approach in the energy range relevant to 
PIXE experimental practice (approximately up to a few tens of 
MeV); it accounts for the energy loss and Coulomb deflection 
of the projectile, the perturbed stationary state and relativistic 
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nature of the target's inner shell. The ECPSSR cross section 
for a given shell is expressed in terms of the PWBA one: 

(rriR ^0 T] \ 
— (C^F — '"'^j 

where C§ is the Coulomb deflection correction, C is the cor- 
rection factor for binding energy and polarization effects, m,R 
is the relativistic correction, go is the minimum momentum 
transfer and 

^ = ..| (6) 

vi being the projectile velocity. 

Further refinements and modifications of the ECPSSR the- 
ory have been proposed: they involve using relativistic Dirac- 
Hartree-Slater wave functions in computing the form factors 
of the theory [24] (identified in the following as ECPSSR-HS), 
account for changes in binding energy of electron states due to 
the presence of the positive ion projectile in the atom (united 
atom correction, identified as ECPSSR-UA), and extend the 
applicability of the theory to relativistic light projectiles for K 
shell calculations [25] (identified as ECPSSR-HE). 

III. Cross section generators 

The evaluation concerns three publicly accessible generators 
of proton ionization cross sections. Two of them, ISICS 
and ERCS08, are distributed through the Computer Physics 
Communications (CPC) Program Library; the third is available 
directly from its author, as specified in [21]. 

A. ISICS 

The ISICS code [18], [26], [27] calculates K, L and M 
shell ionization cross sections by proton and a particle impact 
according to the PWBA approximation and the ECPSSR 
theory in multiple variants. 

The first version of ISICS [18] provided the options of 
PWBA and plain ECPSSR calculations; later versions have 
added capabilites for Hartree-Slater calculations [26], united 
atom option [23] and treatment of relativistic light projectiles 
[27]. The latest version distributed through the CPC Program 
Library at the time of writing this paper is ISICS 2011 [28]; it 
implements a few changes with respect to the previous 2008 
version, that contribute to the numerical robustness of the code. 

ISICS versions up to 2008 used Bearden and Burr's com- 
pilation [29] of atomic binding energies; ISICS 2011 version 
provides also the option of using Williams' compilation [30], 
[31], although the default configuration still uses Bearden and 
Burr's values. 

The theoretical cross sections included in the PIXE data 
library associated with Geant4 9.4 have been produced with 
ISICS 2008. 

B. ERCS08 

ERCS08 [20] is a FORTRAN program for computing 
electron removal cross section by protons and heavier ions. It 
encompasses the calculation of direct ionization and electron 



capture cross sections; only the former is considered in this 
paper. 

The calculations are based on the ECPSSR theory; the 
program allows the configuration of individual components 
of the theory (E, C, PSS, R), as well as other options: R-left 
(as described in [25]), united atom and hSR (as described in 
[24]). 

The program provides a default configuration for each shell 
type (K, L or M), which can be overridden by the user. The 
default configuration for the K shell activates the calculation 
of the form factor involved in the ECPSSR theory using exact 
limits of integration and the correction for relativistic Hartree- 
Slater wave functions for the K shell (the hSR option). The 
united atom approach is taken into account in the default 
configuration for all shells. 

A few test cases documented in |20] show that ERCS08 
calculates ECPSSR cross sections compatible with ISICS 2006 
version to one unit at the fourth significant figure, when it is 
configured to reproduce the same plain ECPSSR settings as 
ISICS 2006. Although [20] states that differences with respect 
to ISICS may be larger when ERCS08 is run using its own 
default parameters and input data, the authors of this paper 
could not find evidence of the experimental vaUdation of the 
default configuration of this code in the Uterature. 

The code includes a set of atomic binding energies; ac- 
cording to [20], they derive from experimental values in the 
literature when available, otherwise from theoretical values 
or, when neither was deemed available, from a first order 
linear regression analysis appUed to available data for elements 
with atomic number greater than 79 documented in [20]. 
Some of the binding energies distributed with ERCS08 could 
be identified as taken from Bearden and Burr's compilation 
[29], the theoretical values of Table IV in [32] and the 1978 
edition of the Table of Isotopes [33]; however, the source 
of other values remains unidentified. Although theoretical 
binding energies are available for any elements with atomic 
number up to 100 in EADL (Evaluated Atomic Data Library) 
[34], they do not appear to have been used in the program. 

C. KIO-LIO 

A software system to calculate ECPSSR cross sections has 
been developed in Pascal by Z. Smit [21] and is available 
on request directly from its author; it can be executed on 
Windows systems. Since it does not appear identified in the 
literature by a specific name, for convenience it is referred to 
in the following as KIO-LIO, as these are the names of the 
executable files provided by the author respectively for K and 
L shell computations. 

The code implements the method described in [21] for the 
K shell and Li, subshells; the cross section for subshell 
ionization is calculated by the same function as that for the 
1/2 subshell, with appropriate binding energy value and further 
multiplied by a factor two. 

The code includes three options of atomic electron binding 
energies, corresponding to the compilations by Bearden and 
Burr [29], the 1978 edition of the Table of Isotopes [33] 
and Sevier [35]; nevertheless, they cannot be selected through 
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the user interface, therefore the distributed executable files 
limit the production to the default configuration with Sevier's 
binding energies. 

D. Comparative features of calculated cross sections 

All of the three generators implement calculations based 
on the ECPSSR theory; nevertheless, as it can be observed 
in Fig. [nil the cross sections they calculate exhibit some 
differences. The figures adopt consistent color, line and marker 
style throughout the paper to facilitate the identification of the 
behavior of each generator 

The histograms in Fig. [T] concern cross sections generated 
by ERCS08 and KIO in their default configurations for K 
shell, and by ISICS configured according to the same nominal 
options as the two other codes: to calculate ECPSSR-HS-UA 
cross sections (nominally equivalent to ERCS08 default) and 
plain ECPSSR ones (nominally equivalent to KIO). The cross 
sections are calculated at the same energies as the experimental 
data collected in [36], which, apart from a few measurements 
at 160 MeV, span the energy range up to approximately 50 
MeV. 

Fig. |2]|4] concerns L shell cross sections calculated at the 
same energies (up to 4 MeV) as the experimental data in [37], 
[38]: by ERCS08 and LIO in their default configurations for 
L shell, and by ISICS configured to calculate ECPSSR-UA 
cross sections (nominally equivalent to ERCS08 default for 
L subshells). The histograms of ISICS plain ECPSSR cross 
sections, corresponding to LIO's default configuration, are not 
superimposed to the other plots for better clarity, since they 
look very close to those produced by the ECPSSR-UA option. 

For nominally equivalent theoretical approaches, discrepant 
results may arise from different mathematical methods adopted 
in the calculations, the use of different atomic parameters (e.g. 
binding energies), different software algorithms or numerical 
precision enforced in the code, or other computational details. 
Comparisons with experimental data are required to quantify 
how the differences of the various cross section generators 
affect the accuracy of the results they produce. 
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Fig. 1. K shell ionization: relative difference of cross sections calculated 
by default ERCS08 configuration (solid black histogram) and KIO (dashed 
red histogram) with respect to ISICS ECPSSR-HS-UA configuration; relative 
difference of cross sections calculated by default KIO configuration with 
respect to ISICS plain ECPSSR configuration (dot-dashed blue histogram). 
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Fig. 2. Li subshell ionization: relative difference of cross sections calculated 
by default ERCS08 configuration (solid black histogram) and LIO (dashed red 
histogram) with respect to ISICS ECPSSR-UA configuration. 



IV. Validation method 

The validation process involves the comparison of cross 
sections calculated by the three generators with experimental 
data. 

A. Experimental references 

The same sources of experimental data used in the validation 
of the Geant4 PIXE data library documented in [14] were 
exploited in the validation process of the three cross section 
generators considered in this paper They are established 
compilations of experimental data in the field; moreover, the 
use of the same experimental references allows a comparative 
evaluation of the accuracy of the newly available generators 
and the currently distributed Geant4 PIXE data library, thus 
facilitating the identification of possible improvements to it. 

The reference experimental data for K shell ionization 
were extracted from the compilation in [36], which includes 
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Fig. 3. L2 subshell ionization: relative difference of cross sections calculated 
by default ERCS08 configuration (solid black histogram) and LIO (dashed red 
histogram) with respect to ISICS ECPSSR-UA configuration. 
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Fig. 4. L-j subshell ionization: relative difference of cross sections calculated 
by default ERCS08 configuration (solid black histogram) and LIO (dashed red 
histogram) with respect to ISICS ECPSSR-UA configuration. 

ionization cross sections along with X-ray and Auger electron 
production cross sections. The production cross sections were 
rescaled into ionization ones using the fluorescence yields 
compiled in [39]; uncertainties were propagated accordingly. 

The reference data for L shell cross sections validation 
were assembled from two complementary collections [37], 
[38]. The experimental compilations report only total L shell 
cross section data for a few lighter elements; sub-shell data 
are listed only for elements with atomic number larger than 
44. The experimental data often exhibit large discrepancies, 
and systematic effects are likely to be present for some 
elements, for which measurements deriving from different 
sources appeal^ to be discordant. 

B. Theoretical cross section production 

ISICS cross sections were calculated using ISICSoo, a 
refactored version of the original ISICS Windows code, which 
has been especially tailored for large scale productions of 
data libraries on Linux platforms. ISICSoo provides the same 
physics functionality as ISICS 2011; its detailed features 
and verification are documented in [40]. The source code 
is intended to be available through the Computer Physics 
Communications Program Library. 

The cross sections generated by ISICSoo appear equivalent 
to those computed by ISICS 2008; the relative differences, 
as shown for example in Fig. |5] for K shell ionization cross 
sections corresponding to plain ECPSSR configuration, are of 
the order of 10^^. 

ERCS08 was compiled and run on a Linux system to pro- 
duce the cross sections subject to validation. The production 
on a Linux platform required a small modification to the 
code released through the CPC Program library, limited to 
the user interface for run control. Input files were generated 
for each of the test cases through the Windows GUI to 
ensure consistency with the default settings provided by the 
system; they activated only the calculation of direct ionization. 
The ERCS08 production system in a Linux environment was 
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Fig. 5. Relative difference of K shell ionization cross sections calculated 
by ISICS 2008 and 2011 version with the plain ECPSSR option; the latter 
were produced by the refactored version of the code identified as ISICSoo. 
The cross sections were computed in both cases at incident proton energies 
and for target elements corresponding to the experimental data in [36]. 

verified to generate identical results to the original code run 
through the Windows GUI, consistent with those documented 
in [20]. 

The production of ISICS and ERCS08 cross sections was 
performed on an Intel Core2 Quad Q9300 CPU computer 
equipped with Scientific Linux 5 operating system and gcc 
4.1.2 compiler 

The executable provided by the author was used for the pro- 
duction of KIO-LIO cross sections on a Microsoft Windows 
XP platform. 

For the purpose of validation the theoretical cross sections 
were calculated at the same energy as the experimental data 
to avoid any perturbation to the intrinsic accuracy of the 
theoretical generators due to interpolation algorithms. This test 
configuration differs from the one reported in [14], where the 
cross sections subject to test were interpolated from tabulations 
produced according to the pre-defined energy grid adopted 
for the data library associated with the Geant4 toolkit. The 
different configuration is justified by the different purpose of 
the test process: validating the Geant4 cross section simulation 
models, that include the interpolation of a data library, in the 
previous paper, and evaluating the intrinsic accuracy of three 
theoretical generators in this one. 

C. Data analysis 

As remarked in the previous sections, the differences among 
the values calculated by the three generators are small (approx- 
imately a few percent); their effect on the accuracy of the cross 
sections can be appreciated only through statistical methods. 

The analysis for the validation of theoretical cross section 
calculations follows the procedure reported in [14]. It is 
articulated over two stages: first a series of tests [41], 
comparing theoretical and experimental cross sections for 
each element and shell or subshell, followed by categorical 
analysis to identify differences across the results of the 
test associated with the three generators. Some experimental 
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Fig. 6. K shell ionization cross sections by proton impact on aluminium 
calculated by theoretical generators and experimental data from [36] (black 
filled circles): ISICS 2011 in ECPSSR-HS-UA configuration (empty circles), 
ERCS08 (empty squares) and KIO (empty triangles), both in default config- 
uration. 



Fig. 7. Li subshell ionization cross sections by proton impact on tellurium 
calculated by theoretical generators and experimental data from [37], [38] 
(black filled circles): ISICS 2011 in ECPSSR-UA configuration (empty 
circles), ERCS08 (empty squares) and LIO (empty triangles), both in default 
configuration. 



data exhibiting large discrepancies with respect to other data 
at the same or neighboring energies, suggesting the presence 
of systematic effects, were excluded from the computation of 
the statistic. The Statistical Toolkit [42], [43] was used for 
goodness-of-fit tests. 

The null hypothesis in the goodness-of-fit testing process 
is defined as the compatibility that the theoretical and experi- 
mental cross section distributions subject to comparison derive 
from the same parent distribution. Unless differently specified, 
a 0.05 significance level is set to define the critical region of 
rejection of the null hypothesis. 

Contingency tables are built from the outcome of the 
test to determine the equivalent behavior of the generators. 
The input to contingency tables derives from the results of 
the test for each element and shell or subshell: they are 
classified respectively as "pass" or "fail" according to whether 
the corresponding p-value is consistent with the defined sig- 
nificance level. The null hypothesis consists in assuming, 
for each contingency table, the equivalence of the generators 
it compares at reproducing experimental measurements. The 
contingency tables are analyzed with Fisher's exact test [44], 
with the test applying Yates' continuity correction [45], 
and with Pearson's test [46] when the number of entries in 
the cells justifies its applicability. A significance level of 0.05 
is set to determine the rejection of the null hypothesis, unless 
specified differently. 

V. Results 

A selection of experimental data and cross sections cal- 
culated by the three generators is displayed in Fig. |6]|9l the 
theoretical values correspond to the generators' configurations 
analyzed in the following sections. 

A comprehensive overview of how theoretical cross sections 
compare to measurements is shown in Fig. [TolfT3] the plots 
display the difference between the cross sections calculated 
by the three generators and experimental data, divided by 
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Fig. 8. L2 subshell ionization cross sections by proton impact on tungsten 
calculated by theoretical generators and experimental data from [37], [38] 
(black filled circles): ISICS 2011 in ECPSSR-UA configuration (empty 
circles), ERCS08 (empty squares) and LIO (empty triangles), both in default 
configuration. 



the corresponding experimental uncertainties. In other words, 
they show by how many standard deviations experimental 
values differ from theoretical ones. These plots include all the 
experimental data of [36]-[38], without excluding any outliers. 
Theoretical generators appear to systematically underestimate 
Li and L2 experimental cross sections. 

The following sections report the statistical analysis to 
evaluate in detail the compatibility of the theoretical cross 
sections with experimental measurements. It is worthwhile to 
note that the incompatibility with experiment of all generators 
in some test cases hints to either an intrinsic deficiency of the 
underlying theory or to systematic effects in the experimental 
data. 
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Fig. 9. L3 subshell ionization cross sections by proton impact on gold 
calculated by theoretical generators and experimental data from [37], [38] 
(black filled circles): ISICS 2011 in ECPSSR-UA configuration (empty 
circles), ERCS08 (empty squares) and LIO (empty triangles), both in default 
configuration. 



Fig. 11. Difference between Li subshell cross sections calculated by the three 
generators and experimental data, divided by the related experimental uncer- 
tainties: ISICS 2011 in ECPSSR-UA configuration (solid black histogram), 
ERCS08 (red dashed histogram) and LIO (dotted blue histogram) in default 
configuration. 
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Fig. 10. Difference between K shell cross sections calculated by the 
three generators and experimental data, divided by the related experimental 
uncertainties: ISICS 2011 in ECPSSR-HS-UA configuration (solid black 
histogram), ERCS08 (red dashed histogram) and KIO (dotted blue histogram) 
in default configuration. 



Fig. 12. Difference between L2 subshell cross sections calculated by the three 
generators and experimental data, divided by the related experimental uncer- 
tainties: ISICS 2011 in ECPSSR-UA configuration (solid black histogram), 
ERCS08 (red dashed histogram) and LIO (dotted blue histogram) in default 
configuration. 



A. K shell 

The p-values of the test comparing calculated and 
experimental K cross sections are reported for each tested 
element in Table U The results are listed for ERCS08 and 
KIO in their default configuration, for ISICS 2011 in various 
options, two of which, ECPSSR-HS-UA and plain ECPSSR, 
are similar to the default configurations of the other codes. 

The results are summarized in Table HH where the test cases 
for which the null hypothesis of compatibility of calculated 
and experimental distributions is rejected at 0.05 significance 
level are counted as "fail", otherwise as "pass". The efficiency 
of each theoretical generator's configuration is defined as the 
ratio between the "passed" and total test cases. 

It is evident from the results in Table HI] that ISICS 2011 
ECPSSR-HS, ECPSSR-HS-UA and ERCS08 default configu- 



ration are the most efficient at reproducing experimental cross 
sections, while KIO and the other ISICS options produce 
slightly less accurate results. Nevertheless, as shown in the 
contingency tables in Table |llll the differences in efficiency 
are not statistically significant at 95% confidence level. 

Although the statistical analysis over the complete data 
samples provides an overall picture of the accuracy of the 
three cross section generators, the detailed results of the 
test in Table U can be useful to optimize the choice in some 
specialized use cases: for instance, ERCS08 cross sections 
for silicon are compatible with experimental data [36], while 
the two other generators fail at reproducing compatible cross 
sections at 0.05 significance level for this element; KIO 
succeeds at reproducing equivalent tungsten cross sections at 
0.05 significance level, while both ISICS ECPSSR-HS-UA and 
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TABLE IV 

P- VALUES OF THE TEST COMPARING CALCULATED AND 
EXPERIMENTAL Li SUBSHELL IONIZATION CROSS SECTIONS BY PROTON 

IMPACT 
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Fig. 13. Difference between Lg subshell cross sections calculated by the three 
generators and experimental data, divided by the related experimental uncer- 
tainties: ISICS 2011 in ECPSSR-UA configuration (solid black histogram), 
ERCS08 (red dashed histogram) and LIO (dotted blue histogram) in default 
configuration. 



ERCS08 default configurations fail. In this respect, however, 
it is worthwhile to note that in the case of tungsten, as well as 
for other test cases concerning heavy elements and relatively 
high 5 — 7 MeV) energy, the Hartree-Slater correction does 
not improve the accuracy of ECPSSR calculations, actually the 
opposite. 

TABLE 111 

Contingency table to estimate the equivalent accuracy of 
ECPSSR K shell cross sections using EADL and Bearden and 
Burr's binding energies 



X"^ test outcome 


ISICS 2011 


ERCS08 


Pass 


51 




51 


Fail 


15 




15 


p-value Fisher test 


1 






p-value Pearson x^ 


1 






p-value Yates x^ 


0.835 




X'^ test outcome 


ISICS 2011 


KIO 


Pass 


51 




47 


Fail 


15 




19 


p-value Fisher test 


0.551 




p-value Pearson x^ 


0.426 




p-value Yates x^ 


0.550 





B. L shell 

The p-values of the test comparing calculated and 
experimental Li, L2 and L3 subshell cross sections are re- 
ported for each tested element in Tables IIVIIVII Whenever the 
data for a given element were present in both experimental 
compilations [37], [38], only the results associated with the 
compilation exhibiting the larger p-values are listed in the 
tables. The results are listed for ERCS08 and LIO in their 
default configuration, for ISICS 2011 in the ECPSRR-UA and 
plain ECPSSR configurations, which are similar to the default 
configurations of the other codes. 

The results are summarized in Table IVIII where the entries 
for "pass", "fail" and efficiency correspond to the same 





ISICS 2011 


ERCS08 


LIO 


z 


ECPSSR 


ECPSSR-UA 


Default 


Default 


45 


0.030 


0.004 


0.021 


0.079 


46 


0.598 


0.528 


0.678 


0.686 


47 


0.966 


0.951 


0.971 


0.979 


48 


0.850 


0.858 


0.865 


0.862 


49 


0.067 


0.087 


0.056 


0.056 


50 


0.197 


0.122 


0.275 


0.344 


51 


0.946 


0.972 


0.936 


0.910 


52 


0.474 


0.419 


0.563 


0.533 


59 


0.056 


0.101 


0.040 


0.022 


62 


0.981 


0.977 


0.992 


0.988 


64 


0.263 


0.243 


0.236 


0.212 


65 


0.974 


0.990 


0.924 


0.865 


66 


0.951 


0.953 


0.939 


0.938 


67 


0.126 


0.127 


0.085 


0.084 


70 


0.619 


0.619 


0.553 


0.552 


72 


0.333 


0.333 


0.418 


0.389 


73 


0.929 


0.929 


0.928 


0.911 


74 


< 0.001 


< 0.001 


< 0.001 


< 0.001 
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0.975 


0.975 


0.986 


0.986 
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0.020 


0.020 


0.003 
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77 


< 0.001 
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0.082 


90 


0.003 


0.003 


0.006 


0.007 


92 


< 0.001 


< 0.001 


< 0.001 


< 0.001 



definitions as in the previous section. ISICS 2011 ECPSSR- 
UA is the most efficient at reproducing experimental cross 
sections, nevertheless, as shown in the contingency tables in 
Table IVIIII the differences in efficiency are not statistically 
significant at 95% confidence level. 

C. Effect of atomic binding energies 

A detailed study [47] devoted to various atomic electron 
binding energy compilations in the context of Monte Carlo 
particle transport highlighted the effect of these atomic pa- 
rameters in the calculation of proton ionization cross sections, 
regarding its accuracy with respect to experimental data. The 
effects are especially visible in the resulting values of K shell 
cross sections. 

As commented in section Ull-Bl ERCS08 includes an ad hoc 
collection of binding energies, that does not correspond to any 
of the compilations used by major Monte Carlo systems and 
specialized PIXE codes. A test was performed to estimate the 
impact of this collection on the accuracy of the cross sections 
by running ISICS with ERCS08 binding energies: a series of 
test compared the compatibility of plain ECPSSR K shell 
cross sections calculated by ISICS based on either its default 
(Bearden and Burr [29]) or ERCS08 binding energies with 
experimental data at 0.05 significance level. The efficiency 
at reproducing experimental data raised from 0.67 ± 0.06 
with Bearden and Burr binding energies to 0.80 ± 0.05 with 
ERCS08 binding energies. 

More extensive documentation of the sources of ERCS08 
binding energies would be beneficial to optimize the accuracy 
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P-VALUES OF THE TEST COMPARING CALCULATED AND EXPERIMENTAL K SHELL IONIZATION CROSS SECTIONS BY PROTON IMPACT 
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13 
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14 
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0.029 
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15 
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0.710 
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0.5 15 


0.695 
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24 


0.374 
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0.1 19 
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0.087 


Jo 


0.019 


0.023 
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0.023 
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50 
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51 


0.787 
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0.787 
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52 
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CI 
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69 


0.014 
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0.015 
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< 0.001 


< 0.001 


< 0.001 


72 


0.275 


0.647 


0.275 
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0.647 


0.650 


0.645 


73 


0.029 


0.004 


0.029 


0.019 


0.004 


0.005 


0.018 


74 


0.232 


< 0.001 


0.232 


0.350 
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0.719 


75 
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0.070 
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0.093 
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0.528 
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83 
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90 
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0.007 


0.003 
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0.001 
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0.646 


0.136 


0.158 


0.646 
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TABLE II 

Summary of the test results of K shell ionization cross sections by proton impact 









ISICS 2011 






ERCS08 


KIO 




ECPSSR 


ECPSSR-HS 


ECPSSR-UA 


ECPSSR-HE 


ECPSSR-HS-UA 


Default 


Default 


Tested elements 


66 


66 


66 


66 


66 


66 


66 


Pass 


44 


51 


44 


46 


51 


51 


47 


Fail 


22 


15 


22 


20 


15 


15 


19 


Efficiency 


0.67±0.06 


0.77±0.05 


0.67±0.06 


0.70±0.06 


0.77±0.05 


0.77±0.05 


0.71 ±0.06 



TABLE Vn 

Summary of test results of L subshell ionization cross sections by proton impact 





ISICS 2011 

ECPSSR ECPSSR-UA 


ERCS08 

Default 


LIO 

Default 




Elements 

Pass 

Fail 

Efficiency 


28 
19 
9 

0.53±0.09 


28 
19 
9 

0.53±0.09 


28 
18 
10 

0.48±0.09 


28 
19 
9 

0.50±0.09 


L2 


Elements 

Pass 

Fail 

Efficiency 


28 
19 
9 

0.68±0.09 


28 
22 
6 

0.79±0.08 


28 
20 
8 

0.71±0.09 


28 
18 
10 

0.64±0.09 


Li 


Elements 

Pass 

Fail 

Efficiency 


28 
25 
3 

0.89±0.06 


28 
25 
3 

0.89±0.06 


28 
26 
2 

0.93±0.05 


28 
21 
7 

0.75±0.08 


L 


Elements 

Pass 

Fail 

Efficiency 


84 
63 
21 

0.75±0.05 


84 
66 
18 

0.79±0.04 


84 
64 
20 

0.76±0.59 


84 
58 
26 

0.69±0.05 



TABLE V 

P-VALUES OF THE X^ TEST COMPARING CALCULATED AND 
experimental L2 SUBSHELL IONIZATION CROSS SECTIONS BY PROTON 



IMPACT 



z 


ISICS 2011 

ECPSSR ECPSSR-UA 


ERCS08 

Default 


LIO 

Default 


45 


0.001 


0.057 


0.023 


< 0.001 


46 


< 0.001 


0.008 


0.002 


< 0.001 


47 


0.015 


0.139 


0.067 


0.029 


48 


0.101 


0.121 


0.080 


0.084 


49 


0.082 


0.256 


0.136 


0.055 


50 


0.092 


0.163 


0.132 


0.073 


51 


1.000 


1.000 


1.000 


1.000 


52 


0.718 


0.722 


0.822 


0.761 


59 


0.350 


0.654 


0.581 


0.258 


62 


0.005 


0.060 


0.026 


0.002 


64 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


65 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


66 


0.836 


0.862 


0.821 


0.777 


67 


0.173 


0.312 


0.196 


0.093 


70 


0.129 


0.139 


0.065 


0.064 


72 


0.984 


0.991 


0.967 


0.962 


73 


0.054 


0.116 


0.065 


0.028 


74 


0.341 


0.453 


0.231 


0.185 


75 


0.398 


0.456 


0.294 


0.283 


76 


0.762 


0.871 


0.677 


0.556 


77 


0.331 


0.433 


0.208 


0.166 


78 


0.655 


0.713 


0.500 


0.478 


79 


1.000 


1.000 


0.998 


0.999 


81 


0.018 


0.018 


0.009 


0.010 


82 


0.841 


0.843 


0.711 


0.728 


83 


0.322 


0.322 


0.187 


0.187 


90 


0.049 


0.049 


0.023 


0.025 


92 


< 0.001 


< 0.001 


< 0.001 


< 0.001 



TABLE VI 

P-VALUES OF THE X^ TEST COMPARING CALCULATED AND 
EXPERIMENTAL L3 SUBSHELL IONIZATION CROSS SECTIONS BY PROTON 



IMPACT 



z 


ISICS 2011 

ECPSSR ECPSSR-UA 


ERCS08 

Default 


LIO 

Default 


45 


0.235 


0.987 


0.948 


0.168 


46 


0.061 


0.829 


0.620 


0.038 


47 


0.436 


0.969 


0.888 


0.347 


48 


0.909 


0.949 


0.919 


0.891 


49 


0.679 


0.978 


0.907 


0.583 


50 


0.221 


0.137 


0.172 


0.204 


51 


0.999 


0.999 


1.000 


0.999 


52 


0.535 


0.463 


0.608 


0.594 


59 


0.019 


0.399 


0.245 


0.010 


62 


0.087 


0.764 


0.586 


0.044 


64 


0.783 


0.274 


0.581 


0.679 


65 


0.309 


0.912 


0.833 


0.188 


66 


0.643 


0.502 


0.589 


0.592 


67 


0.102 


0.512 


0.402 


0.048 


70 


0.599 


0.991 


0.963 


0.394 


72 


0.970 


0.950 


0.992 


0.984 


73 


0.525 


0.867 


0.652 


0.339 


74 


0.701 


0.928 


0.807 


0.549 


75 


0.960 


0.960 


0.992 


0.971 


76 


< 0.001 


< 0.001 


0.004 


< 0.001 


77 


1.000 


0.979 


0.999 


1.000 


78 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


79 


1.000 


1.000 


1.000 


1.000 


81 


0.060 


0.109 


0.066 


0.038 


82 


1.000 


1.000 


1.000 


1.000 


83 


0.460 


0.508 


0.560 


0.430 


90 


0.963 


0.994 


0.978 


0.907 


92 


0.085 


0.021 


0.080 


0.179 
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TABLE VIII 

Contingency table to estimate the equivalent accuracy of 

ECPSSR L SHELL CROSS SECTIONS USING EADL AND BEARDEN AND 

Burr's binding energies 



test outcome 


ISICS 2011 


ERCS08 


Pass 


66 




64 


Fail 


18 




20 


p-value Fisher test 


0.854 




p-value Pearson 


0.712 




p-value Yates 


0.854 




X'^ test outcome 


ISICS 2011 


LIO 


Pass 


66 




58 


Fail 


18 




26 


p-value Fisher test 


0.161 




p-value Pearson x^ 


0.219 




p-value Yates x^ 


0.160 
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of ECPSSR cross section calculations. 

D. Empirical scaling 

Although successful at describing ionization cross sections 
over a wide energy range relevant to PIXE experimental 
applications, the ECPSSR theory fails at reproducing exper- 
imental measurements in some conditions, for instance at 
low energies (below approximately 1 MeV). The Hartree- 
Slater correction applied to the calculation K shell cross 
sections in part overcomes these deficiencies; nevertheless, as 
previously pointed out, this modeling approach exhibits other 
shortcomings, namely for heavy elements at higher energies 
(approximately above a few MeV). 

An empirical correction has been developed by Paul and 
Muhr [48] as a scaling function sj, to be applied to ECPSSR 
cross sections: 

' (7) 



^scaled — 



cecpssr 



The parameters of the scaling function have been fitted to 
experimental data. 

Tabulations of empirically scaled cross sections for K shell 
have been published by Paul and Sacher [36]; the procedure 
of their calculation was modified with respect to that adopted 
in [48], but the updated scaling function is not documented in 
[36]. These tabulations are included in the PIXE data library 
[19]. 

Empirically scaled K shell cross sections are calculated by 
the KIOKC executable in Smit's software system; this program 
uses the atomic binding energies of the 1978 edition of the 
Table of Isotopes, consistently with those used in [48]. 

An implementation of Paul and Muhr's scaling function has 
been developed for use with Geant4 and is intended for release 
in a forthcoming Geant4 version. Using the electron binding 
energies of the 1978 Table of Isotopes, it has been verified 
to reproduce the values reported in [48] for a selection of 
elements and energies. This scaling function can be applied 
to ECPSSR cross sections calculated by other generators to 
improve the calculation accuracy as an alternative to other 
theoretical variants, like the Hartree-Slater and United Atom 
corrections. 

The use of this scaling function could be source of sys- 
tematic effects. The parameters defining the empirical scaling 



Fig. 14. K shell ionization cross sections by proton impact on lead: 
experimental data from [36] (black filled circles), theoretical calculations by 
ISICS 2011 configured with the ECPSSR-HS-UA (empty blues squares) and 
plain ECPSSR (empty black circles) options, and ECPSSR values scaled by 
the empirical scaling function defined in [48] (empty red tiiangles). 
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o 
o 



Energy (MeV) 

Fig. 15. K shell ionization cross sections by proton impact on holmium: 
experimental data from [36] (black filled circles), theoretical calculations by 
ISICS 2011 configured with the ECPSSR-HS-UA (empty blue squares) and 
plain ECPSSR (empty black circles) options, and ECPSSR values scaled by 
the empirical scaling function defined in [48] (empty red tiiangles). 



function documented in [48] derive from a fit to the experi- 
mental data involving calculations of ECPSSR cross sections; 
therefore, strictly speaking, that empirical correction should be 
applied only to identically calculated ECPSSR cross sections. 
The values calculated by the three generators analyzed in 
this paper slightly differ from the ECPSSR cross sections 
reported in [48] for selected elements and energies; however, 
since the discrepancies are relatively small (e.g. less than 1% 
for ECPSSR cross cross sections calculated by ISICS), the 
systematic error due to the application of the empirical scaling 
functions to other ECPSSR calculations than those used in 
[48] would also be small, and negligible with respect to the 
experimental uncertainties of K shell cross sections. 

The contribution of the empirical scaling function to the 
cross section accuracy has been estimated according to the 
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TABLE IX 

Efficiency of empirically scaled and theoretical K shell ionization cross section calculations 







ISICS 2011 




Smit 




ECPSSR 


ECPSSR-HS-UA 


Scaled ECPSSR 


KIO 


KIOKC 


Tested elements 


66 


66 


66 


66 


66 


Pass 


44 


51 


54 


47 


54 


Fail 


22 


15 


12 


19 


12 


Efficiency 


0.67±().06 


().77±0.05 


0.82±0.05 


0.71±0.06 


0.82±().05 



procedure described in section IIV-CI this process does not 
constitute a validation, since the experimental data to which 
the calculated cross sections are compared are in large part the 
same used for fitting the parameters of the scaling function 
itself. The usual procedure to account for fitting constraints 
in a -y^ test, consisting of reducing the number of degrees 
of freedom accordingly, cannot be applied to this case due 
to incomplete knowledge of the fit configuration. In fact, the 
fits described in [48] concern groups of several elements, 
while in this analysis the tests are performed for individual 
elements; moreover, the experimental samples of Paul and 
Sacher's paper [36], used in this analysis, do not exactly 
coincide with the data on which the scaling function drawn 
from Paul and Muhr's paper [48] has been fitted, although a 
large fraction are common. The ECPSSR cross sections, which 
are involved in the fit, are slightly different in [48] and in this 
analysis. Therefore the results reported in the following should 
be considered only as a demonstration of the capability of the 
empirical scaling to describe the data. 

The efficiency, defined as in the previous sections, is listed 
in Table UX] for ECPSSR cross sections calculated by KIOKC 
and by ISICS scaled by Paul and Muhr's empirical function 
along with the efficiency of the ISICS ECPSSR and ECPSSR- 
HS-UA options, and of KIO. The gain in efficiency due to the 
scaling function appears substantial, although, as illustrated 
in the contingency table reported in Table |X] the outcome 
of the statistical analysis is controversial: the hypothesis of 
equivalence of theoretical and empirically scaled ECPSSR 
cross sections is rejected with 0.05 significance by Pearson's 
test, but only with 0.1 significance by Yates' test and 
Fisher's test. The gain with respect to the theoretical ECPSSR- 
HS-UA option is not statistically significant; nevertheless, 
in some experimental scenarios for which the theory does 
not reproduce the measurements accurately, for instance as 
illustrated in Fig. [14] the capability of applying empirical 
scaling corrections to the theoretical cross sections could be a 
valuable option. However, within the data sample analyzed in 
this paper, one can also identify cases where theoretical cross 
sections calculated with the Hartree-Slater and United Atom 
corrections are more accurate than empirically scaled ECPSSR 
ones; an example is shown in Fig. [15] 

E. Interpolated cross sections 

In the data driven approach to PIXE simulation described 
in [14], the cross sections calculated by analytical generators 
are tabulated at predefined energies and used in Monte Carlo 
simulation to compute cross sections corresponding to the 
energy of a given particle being tracked through matter. 



TABLE X 

Contingency table to estimate the equivalent accuracy of 
plain and empirically scaled isics 201 1 ecpssr k shell cross 
sections 



test outcome 


ECPSSR 


Scaled 


Pass 


44 


54 


Fail 


22 


12 


p-value Fisher test 


0.072 


p-value Pearson 


0.047 


p-value Yates 


0.073 



A test was performed to verify if the interpolation process 
performed in Geant4 would significantly affect the accuracy 
of the theoretical cross sections. 

The distribution of the relative difference between inter- 
polated and theoretical ECPSSR cross sections is shown in 
Fig. [16] for K shell cross sections; the distributions for L 
subshells are similar. The data span the energy range covered 
by experimental measurements; the theoretical cross sections 
are generated by ISICS, while the interpolated values are com- 
puted by logarithmic interpolation of the tabulations produced 
by ISICS. The plot includes only truly interpolated values: 
cross sections calculated by the interpolation algorithm at 
energies corresponding to the tabulated values are excluded, 
since the differences are obviously null. The distribution can be 
described by a Gaussian with 0. 17% mean and 0. 1 1 % standard 
deviation; the tail corresponds to regions of relatively steep 
variation of the cross sections as a function of energy, mostly 
occurring with light targets at low energies (below 1 MeV). 
The error due to interpolation is negligible with respect to 
the experimental uncertainties of the cross sections, as can be 
observed in Fig. [17] 

Interpolated cross sections exhibit the same compatibility 
with experiment as theoretically calculated ones: the outcome 
of the x^ test, as rejection or non-rejection of the null hypoth- 
esis at 0.05 significance level, is strictly identical for ECPSSR 
cross sections calculated by ISICS and interpolated from 
tabulations produced with ISICS, when they are compared to 
the same experimental samples applying the same analysis 
criteria, although the test statistic assumes slightly different 
numerical values for the two categories in the individual 
test cases. This result holds for cross sections interpolated 
from ISICS 2008 tabulations with respect to theoretical values 
directly calculated by either ISICS 2008 or 2011. 

From this test one evinces that the interpolation process 
used in Geant4 PIXE simulation [14] introduces a negligible 
perturbation to the accuracy of the cross sections calculated in 
the course of the simulation; it also confirms that the updates 
implemented in ISICS 2011 do not modify the accuracy 
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Fig. 16. Absolute relative difference between interpolated and theoretical 
ECPSSR K shell ionization cross sections by proton; the theoretical cross 
sections are produced by ISICS, while the interpolated values are calculated by 
the logarithmic interpolation algorithm encompassed in the data management 
component described in [14], based on tabulations produced by ISICS. 
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Fig. 17. Ratio of K shell cross section experimental uncertainties and 
interpolation errors; the interpolation error is defined as the absolute difference 
between interpolated and theoretical cross sections. The theoretical cross 
sections are produced by ISICS, while the interpolated values are calculated by 
the logarithmic interpolation algorithm encompassed in the data management 
component described in [14], based on tabulations produced by ISICS; the 
experimental uncertainties are taken from [36]. 



of Geant4 PIXE simulation with respect to the data library 
distributed with Geant4 9.4, which is based on the previous 
ISICS version. 

It is worthwhile to stress that these results concerning 
the interpolation of cross section tabulations produced by 
ISICS concern the data library and interpolation algorithm 
described in [14]; they would not necessarily hold for tabula- 
tions adopting different sampling intervals or exploiting other 
interpolation algorithms. 

VI. Conclusion 

Three pubUcly available theoretical generators of inner shell 
ionization by proton impact have been analyzed to assess the 
state-of-the-art in the field: the 2011 version of ISICS, the 



new code ERCS08 and Smit's software system (here called 
KIO-LIO). 

The accuracy achieved by the three generators has been 
estimated through the statistical comparison of their results 
with experimental K and L shell cross sections. The analysis 
confirms quantitatively the qualitative similarity of the three 
generators. For K shell ionization, ISICS and ERCS08 exhibit 
the same overall efficiency at reproducing experimental ref- 
erences, while the KIO system is sUghtly less accurate. For 
the L shell, ISICS achieves the highest efficiency and LIO the 
lowest: nevertheless, the differences among the generators are 
not statistically significant at 0.05 level. One can conclude that 
ISICS represents the state-of-the-art among freely available 
ECPSSR generators. 

Although quite similar in their overall statistical perfor- 
mance, the three generators exhibit distinctive features for 
a few specific elements (for instance, for siUcon K shell 
ionization). 

Their documentation in the detailed results of this paper 
provides guidance to optimize the selection of cross section 
options in applications, which may be sensitive to the accuracy 
of the results for specific targets of experimental interest. 

These results are relevant to the production of cross sec- 
tion tabulations for Monte Carlo particle transport, and the 
selection of optimal cross section options in simulation appli- 
cations. 

The analysis also assessed the equivalent compatibility with 
measurements of the 2011 and 2008 versions of ISICS; the 
latter was used for the production of the cross section data 
library released in Geant4 9.4 for PIXE simulation. Moreover, 
it demonstrated that the interpolation of tabulated values does 
not significantly affect the accuracy of the cross sections 
calculated in Geant4. 

The PIXE data library has been extended to encompass 
tabulations produced by the other generators analyzed in this 
paper in addition to those based on ISICS, so that users can 
profit from the distinctive features of the three generators 
identified in this study. An updated version is intended to be 
released with a forthcoming version of the Geant4 toolkit and 
through the RSICC distribution center. An extension to Geant4 
code is also foreseen to provide functionality for empirically 
scaling ECPSSR K shell cross sections. 

Apart from the calculation of cross sections for the ion- 
ization of target atoms, which is the subject of this paper, 
the simulation of PIXE with Geant4 involves the atomic 
relaxation, whose accuracy is documented in [49]. Other issues 
concerning the consistent treatment of the discrete process of 
atomic relaxation along with the ionization process, which is 
affected by infrared divergences, are extensively discussed in 
[14]. 
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